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O ' Abstract 

' One way to ameliorate the SUSY flavor and CP problems is to postulate that scalar 

i ' masses lie in the TeV or beyond regime. For example, the focus point (FP) region of 

the minimal supergravity (mSUGRA) model is especially compelling in that heavy scalar 
masses can co-exist with low fine-tuning while yielding the required relic abundance of 
cold dark matter (via a mixed higgsino-bino neutralino). We examine many of the charac- 
teristics of collider events expected to arise at the CERN LHC in models with multi-TeV 
scalars, taking the mSUGRA FP region as a case study. The collider events are charac- 
terized by a hard component arising from gluino pair production, plus a soft component 
arising from direct chargino and neutralino production. Gluino decays in the FP region 
are characterized by lengthy cascades yielding very large jet and lepton multiplicities, and 
a large 6-jet multiplicity. Thus, as one steps to higher jet, 6-jet or lepton multiplicity, 
signal-over-background rates should steadily improve. The lengthy cascade decays make 
mass reconstruction via kinematic edges difficult; however, since the hard component is 
nearly pure gluino pair production, the gluino mass can be extracted to ±8% via total 
rate for E™ lsSj t- > 7-jet + > 2 fo-jet events, assuming 100 fb _1 of integrated luminosity. 
The distribution of invariant mass of opposite-sign/same-flavor dileptons in the hard com- 
ponent exhibits two dilepton mass edges: — m~ and — ■ As a consistency 
check, the same mass edges should be seen in isolated opposite-sign dileptons occurring 
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in the soft component trilepton signal which originates mainly from chargino-neutralino 
production. 
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1 Introduction 



The minimal supergravity (mSUGRA) model[TJ[2j[3] is a well-motivated supersymmetric model 
with a small parameter space that forms a template for many investigations of the phenomeno- 
logical consequences of weak scale supersymmetry. In mSUGRA, it is assumed that the minimal 
supersymmetric standard model, or MSSM, is a valid effective theory of physics between the 
energy scales Q = M GUT and Q = M weak . It is further assumed that all MSSM scalar masses 
unify to a common value mo at Mqut, while gauginos unify to a common value mi/2 and trilin- 
ear soft terms unify to a common value Aq. The weak scale soft parameters can be calculated 
by renormalization group evolution from Mqut to M wea k- The large value of the top quark 
Yukawa coupling drives the up-Higgs squared mass to negative values, leading to radiative elec- 
troweak symmetry breaking (EWSB). The EWSB minimization conditions allow one to trade 
the bilinear soft breaking term B for tan f3, the ratio of Higgs field vevs. It also determines the 
magnitude (but not the sign) of the superpotential Higgs mass term /i. Thus, the entire weak 
scale sparticle mass spectrum and mixings can be calculated from the well-known parameter 
set 

m , mi/2, A , tan/5, signiji). (1) 

Thus, once this parameter set is stipulated, a whole host of observables, including the neutralino 
dark matter relic density Q-g h 2 and collider scattering events, may be calculated. For imple- 
mentation, we use Isajet v7.74[U [5] to calculate the sparticle mass spectrum and associated 
collider events, and IsaReD[6] to calculate the neutralino relic density 

One of the important consequences of the MSSM, due to -R-parity conservation, is that 
the lightest SUSY particle (LSP) is absolutely stable. In mSUGRA, the LSP is usually found 
to be the lightest neutralino Zi, which is a weakly interacting massive particle (WIMP), and 
hence has the potential to naturally match the measured abundance of cold dark matter in the 
universe. An analysis of the three- year WMAP and galaxy survey data sets [7] implies that the 
ratio of dark matter density to critical density, 

Vzh 2 = PzJp c = OAll±™\l (2a). (2) 

where h = 0.74 ± 0.03 is the Hubble constant. By comparing the mSUGRA predicted value 
of h 2 t° this measured value, one finds that only certain parts of the mSUGRA parameter 
space are cosmologically allowed. These include the following. 

• The bulk region at low m and low mm, where neutralino annihilation is enhanced by 
light t-channel slepton exchange [8j. The tight WMAP QcDAdh 2 limit has pushed this 
allowed region to very small mo and m\ji values, while LEP2 limits on and m^ (and 
possibly on m/,) have excluded these same low values so that almost no bulk region has 
survived [9]. 

• The stau co-annihilation region occurs at very low mo but any m^ values, so that 

~ nig , and neutralinos can annihilate against tau sleptons[10j in the early universe. 
For certain A values which dial mi to very low values, there also exists a top-squark 
co-annihilation region [TT]. 
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• The A-funnel region occurs at large values of the parameter tan/3 ~ 50, where 2m^ ~ 
vnAi and neutralinos can annihilate through the broad pseudoscalar Higgs resonance [12] . 
There is also a light Higgs resonance region where 2m^ ~ rrih at low mi/2 values[Hl [13]. 

• At large mo near the boundary of parameter space, the superpotential Higgsino mass 
term \i becomes quite small, and the Z\ can become a mixed higgsino-bino neutralino. 
This region is known as the hyperbolic branch/focus point region (HB/FP) [TU HSl HSj. 
In this case, neutralino annihilation to vector bosons is enhanced, and a match to the 
WMAP measured relic density can be found. 

The HB/FP region of mSUGRA is especially compelling. In this region, the large value 
of ttiq ~ several TeV means that possible SUSY contributions to various flavor- changing and 
CP- violating processes are suppressed by the large squark and slepton masses. For instance, 
SUSY contributions to the flavor- violating decay b — > 57 are small, so in the HB/FP region 
this decay rate is predicted to be in accord with SM predictions, as observed. Meanwhile, the 
calculated amount of fine-tuning in the electroweak sector has been shown to be small [T4"[ [T5] , 
in spite of the presence of multi-TeV top squarks. 

In this paper, we examine the HB/FP region with regard to what sort of collider events are 
expected at the CERN LHC pp collider, which is set to begin operating in the near future at 
a center-of-mass energy y/s = 14 TeV. Much previous work on this issue has been done. In 
Ref. [TTl [TS], the reach of the LHC in the mSUGRA model, including the HB/FP region, was 
calculated. The reach for 100 fb _1 was found to extend to rai/2 ~ 700 GeV, corresponding 
to a reach in m- g of about 1.8 TeV. The reaches of y/s = 0.5 and 1 TeV e + e~ linear colliders 
were also calculated [19], and found to extend past that of the LHC, since when /x becomes 
small, charginos become light, and chargino pair production is a reaction that e + e~ colliders 
are sensitive to, essentially up to the kinematic limit for chargino pair production. In fact, 
the reach of the Fermilab Tevatron for SUSY in the clean trilepton channel [2Tj is somewhat 
enhanced in the HB/FP region[23], since charginos and neutralinos can be quite light, and 
decay with characteristic dilepton mass edges. The reaches of direct [25] and indirect [26] dark 
matter search experiments are also enhanced in the HB/FP region. 

In Ref. [27], Hinchliffe and Paige examined characteristic measurements that the LHC 
could make for an mSUGRA sample point nearby to the HB/FP region. They found a good 
signal/background ratio could be obtained with a hard cut on effective mass M e ff = Yl& ts Et + 
J2ie P tons Et + E™ 1SS (e.g. M e ff > 400 GeV) and by requiring the presence of a 6-jeto Some 
characteristic distributions such as m(£ + £~) which gave a dilepton mass edge at ~ m z 1 anc ^ 

m(b — jet, i) < J (mf — M^)/2 (indicating the presence of a t-quark in the decay chain) could 
be made. 

In Ref. [29], Mizukoshi, Mercadante and Tata found that the LHC reach in the HB/FP 
region could be enhanced by up to 20% by requiring events with the presence of one or two 
tagged fo-jets. In Ref. [30J, a model-independent exploration of the HB/FP region was made 
with regard to collider and dark matter signals. The LHC reach via multilepton cascade decays 
was compared to the LHC reach via clean trileptons from pp —>■ W\Z 2 — > 3£ + E™ lss production. 
In the latter process, backgrounds from W*Z* and W*^* were calculated, and the trilepton 

1 Thc effective mass was introduced and used in heavy top quark production [28] . 
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reach was found to be comparable to- but slightly smaller than- the reach via a search for gluino 
cascade decays. In Ref. [HT], the authors examined what sort of cosmological measurements 
could be made in several mSUGRA case studies (including the pointl§ LCC2 in the HB/FP 
region) by measurements at the LHC and a y/s = 0.5 and 1 TeV ILC. For LCC2 at the LHC, 
they assumed the Z 3 — Z\ and Z 2 — Z\ mass edges could be measured to an accuracy of 1 GeV, 
while it was conjectured that rrig and m^r could be measured to ~ 10% accuracy via some 
kinematic distributions. 

In this paper, in anticipation of the LHC turn on, we wish to understand many of the 
characteristics of collider events expected in the HB/FP region, with an eye towards sparticle 
mass measurements rather than reach studies. We find that the expected collider events in the 
HB/FP region separate themselves into a hard component, arising from gluino pair production, 
and a soft component, arising from pair production of charginos and neutralinos. The gluino 
pair production events typically involve lengthy cascade decays to top and bottom quarks [32] , 
and so high jet, 6-jet and isolated lepton multiplicities are expected. However, the complex 
cascade decays do not lend themselves to simple kinematic measurements of the gluino or 
neutralino masses, mainly due to the combinatorics of picking out the correct gluino decay 
product^). We do find that the gluino mass should be extractable based on total rate in the 
multi-jet+multi-lepton+£'™ lss events to a precision of about 5-10% for 100 fb _1 of integrated 
luminosity. For both the hard and soft components, the Z% — Z\ and Z 2 — Z\ mass edges should 
be visible. 

The remainder of this paper is organized as follows. In Sec. [21 we present details of sparticle 
masses and cross sections expected from the HB/FP region at the LHC. In Sec. El we present 
some details of our signal and background calculations. In Sec. S]we present distributions for 
a variety of collider observables for a case study and SM backgrounds. We present a set of 
cuts that allows good separation of signal vs. background over a large range of nig values. In 
Sec. [51 we show expected signal-to-background plots for gluino pair production and discuss 
how these can be used to extract a measurement of the gluino mass. In Sec. El we address 
leptonic signals. We conclude in Sec. [3 



2 Sparticle production and decay in the HB/FP region 

In the mSUGRA model, for a given set of GUT scale soft SUSY breaking (SSB) masses, the 
associated weak scale values may be computed via renormalization group (RG) evolution |3~5] . 
Once the weak scale SSB terms have been obtained, then the scalar potential must be minimized 
to determine if electroweak symmetry is properly broken. While one EWSB condition allows 
the bilinear parameter B to be traded for tan f3, the other condition reads (at one-loop) 

2 _ m H d -m 2 Hu tan 2 f3 _ M% 
^ (tan^-l) 2 ' 1 } 

which determines the magnitude of the superpotential /1 parameter. Thus, one condition that 
EWSB is successfully broken is that a positive value of /j, 2 has been generated. Roughly, if all 

2 The parameters of the point are m = 3280 GeV, m 1 / 2 = 300 GeV, tan/3 = 10, A = GeV and /i > 0. 
3 Studies of gluino mass |33| and spin |34| determination have been made through the cascade decays g — > bb\ 
with 61 — ► Z2 — > I — > Z\ at parameter point SPSla. 
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the soft parameters entering Eq. [3] are of order M|, then naturalness is satisfied, and the model 
is not fine-tuned. 

For a fixed value of the parameter m\/2 in the mSUGRA model, if mo is taken to be of order 
the weak scale, then m 2 Hu is driven to negative values at the weak scale owing to the push from 
the large top quark Yukawa coupling in the RGEs. However, if is taken too large, then the 
GUT scale value of m 2 H is so high that it is not driven to negative values when the weak scale 
is reached in RG running, and a positive value of /x 2 cannot be found. Intermediate to these 
two extreme cases must exist a region where /j, 2 is found to be zero, which forms the large m 
edge of parameter space. If fi 2 is positive, but tiny, then extremely light higgsino-like charginos 
will be generated, in conflict with bounds from LEP2, which require > 103.5 GeV. If fi 2 
is large enough to evade LEP2 limits, then large higgsino-bino mixing occurs in the chargino 
and neutralino sectors, and in fact the lightest neutralino becomes a mixed higgsino-bino dark 
matter particle. A lightest neutralino of mixed higgsino-bino form has a large annihilation rate 
to vector bosons in the early universe, and hence may have a dark matter relic density in accord 
with WMAP measurements. In this region, dubbed the hyperbolic branch/focus point region, 
multi-TeV squark and slepton masses can co-exist with low fine-tuning as dictated by Eq. [3j 
Thus, the HB/FP region is characterized by TeV-scale squark and slepton masses, which are 
useful for suppressing possible FCNC or CP-violating processes, low fine-tuning, and a dark 
matter relic density in accord with WMAP. Given these qualities, it is important to investigate 
what HB /FP supersymmetry events would look like at the LHC collider and what sort of mass 
measurements could be made in this region. 
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Table 1: Points in the HB/FP region that yield a relic density £1% h 2 ~ 0.11. Common values 
of tan/3 = 30, Aq = 0, fi > and m t = 175 GeV are assumed. The gluino mass in GeV and 
its anticipated measurement uncertainty at the LHC are also given. The total gluino width is 
given in MeV. 

We have generated sparticle mass spectra in the HB/FP region using Isajet v7.74, retaining 
only points which yield a relic density VL~ h 2 ~ 0.11. These points and the anticipated measure- 
ment uncertainty are listed in Table [TJ Since the scalar quarks are decoupled, the gluino decays 
via suppressed 3-body decays, making the gluino width small and not a significant source of 
mass measurement uncertainty. In Fig. [TJ, we show the points in the HB/FP from Table [TJ in 
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Figure 1: Regions of allowed relic density in the m vs. mi/ 2 plane for A = 0, tan/3 = 30 and 
m t = 175 GeV. Points in the HB/FP region from Tabled] that are consistent with the observed 
relic density ~ 0.11 are shown by a blue x. The white region above the focus point and 

surrounding the coannihilation region give and overdensity of neutralino dark matter. 



the mo vs. mi/2 plane. Isajet uses two-loop RGEs for the scalar mass evolution, and minimizes 
the RG-improved one-loop effective potential at an optimized scale Q = y/rrit L m t R (which ac- 
counts for leading two-loop effects). A unique feature of Isajet's sparticle mass algorithm is 
that it decouples various SSB terms from the RG evolution at their own mass scales, which 
gives a more gradual transition from the MSSM to the SM effective theory, as opposed to other 
approaches which use an "all-at-once" transition. Thus, the Isajet algorithm should give a good 
representation of sparticle mass spectra in cases that involve a severely split mass spectrum, 
such as in the HB/FP region. 

In Fig. [21 we show the sparticle mass spectra as a function of m.1/2 along a line with 
^Zi^ 2 ~ 0-11? f° r tan/5 = 30, Aq = and fi > 0. The physics in the HB/FP region is not very 
sensitive to tan/3 or A , since the scalar masses effectively decouple. We take m t = 175 GeV, 
but note that the mo value needed to obtain the correct relic density is extremely sensitive to 
the value of m t used, as shown in Ref. [23]. In our case, since the scalar masses are expected 
to decouple, the m t dependence should not matter greatly for the phenomenology of interest]^ 

While squarks, sleptons and heavy Higgs scalars range in mass from 2.5 — 4.5 TeV along 
the range of mi/2 shown in Fig. [21 the g remains relatively light, of order 650 — 2200 GeV. 
In addition, since /i and are low, the charginos and neutralinos are all quite light, and 
possibly accessible to LHC experiments. The lower edge of the plot where m 1( / 2 ~ 250 GeV is 
excluded by the LEP2 constraint on the chargino mass. 

In Fig. [31 we show sparticle pair production rates as a function of m- g in the HB/FP 
region. While the production cross sections are evaluated at lowest order in perturbation 

4 In our study, we adopt the reference value m t = 175 GeV to allow comparisons with other studies. The 
recent world average for the t-quark mass is m t = 171.4 GeV pH] . 
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HB/FP Sparticle Masses 




theory, we adopt a renormalization/factorization scale choice Q = (mi + 77*2) /4 for the gluino 
pair production cross section which gives good agreement between LO and NLO results [36] 1^1 
For low values of m g ~ 700 GeV, gluino pair production is in the pb range, while a variety 
of chargino and neutralino production processes (e.g. W\Z\ t 2,z and W^Wx production) have 
comparable rates. For higher values of rrig, the gluino pair production cross section drops 
quickly, and is below the fb level for m g > 1900 GeV. The various chargino and neutralino 
production rates drop less quickly, and turn out to be by far the dominant sparticle production 
cross sections for m g ~ 1.5 TeV. 

In the WMAP- allowed HB/FP region, since squarks have masses in the TeV range, only 
three-body decay modes of the gluino are allowed. Moreover, since the \x parameter is small 
and the lighter inos have a large higgsino component, the third generation quark- squark-ino 
couplings are enhanced by top quark Yukawa coupling terms [3"T[ 13"%] . and gluinos dominantly 
decay to third generation particles, especially the top quark. Thus, the dominant gluino decays 
modes in the HB/FP region consist of g — > ttZj or tbWj. Some major g branching fractions are 
listed in Table [2] for a case study which we label as FP5 with m g = 1076 GeV. The Feynman 
diagrams of these dominant decay modes are shown in Fig. HI Thus, we expect in the HB/FP 
region that pp — > ggX will yield events with very large jet and 6-jet multiplicities, and isolated 
leptons. However, the combinatoric backgrounds will likely make kinematic reconstruction 
of mass edges which depend on very difficult. Meanwhile, for the same case study as in 
Table [21 since Z2 — > e r e~Z\ and Z3 — > e + e~Zx both occur at a branching fraction of 3.4%, it 
might be possible to see both the Z 2 — Z\ and Z 3 — Z\ mass edges in distributions of invariant 

5 NLO gluino. chargino and neutralino cross sections are shown versus weak scale gaugino mass M\ in the 
HB/FP region in Ref. [30]. 
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Figure 3: Selected sparticle pair production cross sections vs. m~ g along a line of constant 
fi~ h 2 = 0.11 in the HB/FP region of mSUGRA with A = 0, tan/5 = 30, \i > and m t = 175 
GeV. 




3 HB/FP signal and background event generation 

We use Isajet 7.74(1] for the simulation of signal and background events at the LHC. A toy 
detector simulation is employed with calorimeter cell size At] x A<p = 0.05 x 0.05 and —5 < 
rj < 5. The HCAL energy resolution is taken to be 80%/v^E + 3% for \rj\ < 2.6 and FCAL is 
100^0/^+5% for > 2.6. The ECAL energy resolution is assumed to be 3%/\^E+0.5%. We 
use a UAl-like jet finding algorithm with jet cone size R = 0.4 and require that E-rijet) > 50 
GeV and \rj(jet)\ < 3.0. Leptons are considered isolated if they have pr{e or fi) > 20 GeV and 
\rj\ < 2.5 with visible activity within a cone of AR < 0.2 of E.Ef? < 5 GeV. The strict isolation 
criterion helps reduce multi-lepton backgrounds from heavy quark (cc and bb) production. 
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Table 2: Selected branching fractions of the g for FP5 case study with parameters m = 3050 
GeV, m V 2 = 400 GeV, A = 0, tan/5 = 30 and // > 0. 
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QCD (p T : 200 - 
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4.4 


W + jets; W -> 


e,fi,T (p T (W) 


: 100 - 4000 GeV) 


5 x 10 5 


3.9 x 10 5 


1815.9 


Z + jets; Z — > rr, i^s (pr(Z) : 


100 - 3000 GeV) 


5 x 10 5 


1.4 x 10 5 


845.3 


tt 






3 x 10 6 


4.6 x 10 5 


6415.8 


WW, zz, wz 






5 x 10 5 


8.0 x 10 4 


9.3 


signal (FP5: m- g 


= 1076 GeV) 




2 x 10 5 


1.2 x 10 3 


77.5 



Table 3: Events generated and cross sections for various SM background processes plus one 
HB/FP case study FP5 with m = 3050 GeV, m 1/2 = 400 GeV, A = 0, tan/3 = 30 and [i > 0. 
The CI cuts are specified in Eqns. (4 — 7). 



We identify a hadronic cluster with Et > 50 GeV and \rj(j)\ < 1.5 as a fo-jet if it contains a 
B hadron with Pt(B) > 15 GeV and |?7(-B)| < 3 within a cone of AR < 0.5 about the jet axis. 
We adopt a fo-jet tagging efficiency of 60%, and assume that light quark and gluon jets can be 
mis-tagged as 6-jets with a probability 1/150 for E T < 100 GeV, 1/50 for E T > 250 GeV, with 
a linear interpolation for 100 GeV< Et < 250 GeV. 

We have generated 200K events each for a variety of min values in the HB/FP region 
restricted to have ^^ i h 2 ~ 0.11. In addition, we have generated background events using Isajet 
for QCD jet production (jet-types include g, u, d, s, c and b quarks) over five pr ranges as 
shown in Table [3j Additional jets are generated via parton showering from the initial and 
final state hard scattering subprocesses. We have also generated backgrounds in the W + jets, 
Z + jets, (175) and WW, WZ, ZZ channels at the rates shown in Table [3J The W + jets 
and Z + jets backgrounds use exact matrix elements for one parton emission, but rely on the 
parton shower for subsequent emissions. 
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4 Event characteristics in the HB/FP region 

We begin by applying a set of pre-cuts to our event samples, which we list as cuts set CI [22 



CI Cuts: 



771 miss 

Ulrp 

n(jets) 
ErUl,j2,j3,j4) 
St 



> 



> 



> 



> 



(100 GeV,0.2M e// ), 
4, 

100, 50,50,50 GeV, 
0.2. 



(4) 
(5) 
(6) 
(7) 



Here, M eff is defined as in Hinchliffe et a/.[22] as M eff = E^ 88 + E T (jl) + E T (j2) + E T (j3) + 
_&rC/4), where jl — j4 refer to the four highest Et jets ordered from highest to lowest Et, 
_E™ 1SS is missing transverse energy and St is transverse sphericity. The event rates in fb are 
listed after CI in Table [31 and we find that signal with these cuts is swamped by various SM 
backgrounds, especially those from QCD multi-jet production and ti production. 

Next, we investigate a variety of distributions. We show in Fig. [5] the M e ff distribution 
after using CI. The gray histogram denotes the sum of all backgrounds, while individual BG 
contributions are identified by the legend. The signal for case study FP5 is denoted by the purple 
histogram. In many models investigated by Hinchliffe et ai, it was found that signal emerges 
from BG at an M e ff value near the peak of the distribution, which in fact provides a rough 
estimate of the strongly interacting sparticle masses involved in the production subprocess. In 
the HB/FP region, however, squarks have decoupled from the hadronic sparticle production 
cross section, so only gluino pair production contributes. In addition, since in the HB/FP region 
gluinos decay via three-body modes, the average jet E T is reduced significantly compared to 
SUSY cases with similar sparticle masses but with dominantly 2-body decays. Hence, in the 
HB/FP region, the M e ff distribution from the signal is typically buried beneath SM BG. In 
addition, for this case study, we see some structure to the M e ff distribution in the form of 
two separate peaks (which stand out more clearly on a linear scale, when BG is neglected). 
The peak near M e ff ~ 500 GeV comes dominantly from the soft signal component, which is 
mainly high pt chargino and neutralino production, which after all is the dominant sparticle 
production process in the HB/FP region. A second peak around M e ff ~ 1200 GeV comes from 
gluino pair production, which we denote as the hard component of the signal. 

We noted earlier, based on an examination of gluino decay modes in the HB/FP region, that 
LHC collider events ought to be characterized by large jet multiplicity, large fr-jet multiplicity 
and large isolated lepton multiplicity. With this in mind, we show in Fig. [6] the multiplicity 
of jets expected from signal and from SM BG, after cuts CI. At low n(jets) ~ 4 — 6, the 
distribution is dominated by QCD, ti and W, Z +jets production. However, at much higher jet 
multiplicities ~ 9 — 10, the signal distribution emergeqj from the BG. Of course, at these high 
jet multiplicities, one may question the validity of the theoretical BG calculations. However, 
by investigating QCD multijet production and W, Z + jets production without imposing CI, it 

6 Thc use of the steps in the jet multiplicity was introduced in Ref. |39j in extracting the signal of top quark 
pair production. 
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GeV, A = 0, tan/? = 30, ji > and m t = 175 GeV (where rrig = 1076 GeV), versus various 
SM backgrounds. 



may be possible to normalize the expected BG distributions to measured data, and thus obtain 
after LHC turn-on improved estimates of expected BGs in these channels. 

In Fig. [7J we show the expected multiplicity of fe-jets for signal and SM BG. The soft 
component of signal is expected to be 6-jet poor, since it comes from hadronic chargino and 
neutralino decays. However, the hard component is expected to typically contain at least 4 
6-jets, aside from efficiency corrections. Indeed, we see that the signal distribution extends out 
to high 6-jet multiplicities of n(b — jet) ~ 5 — 8, while the BG typically gives — 2 6-jets. As 
noted previously, Mercadante et al. exploited this fact to enhance the LHC reach for SUSY in 
the HB/FP region [29]. 

In Fig. [HI we show the multiplicity of isolated leptons: electrons or muons. Again, while 
low lepton multiplicity is dominated by SM backgrounds, the high lepton multiplicity should be 
dominated by signal, owing to the lengthy gluino cascade decays, which can spin off additional 
isolated leptons at various stages. 

At this point, it is evident that requiring collider events with high jet and high b-jet multi- 
plicity will aid in separating signal from BG in the HB/FP region. Thus, in Fig. [91 we show 
the augmented effective mass distribution At, where 

A T = E™ S + J2 E t + Y, E t, (8) 

leptons jets 

which gives the added contribution of additional jets beyond n(jets) = 4 and also a contribution 
from isolated leptons. The distributions in Fig. Oall contain, along with cuts CI, n(jets) > 6 
and a) n(b — jets) > 0, b) n(b — jets) > 1, c) n(b — jets) > 2 and d) n{b — jets) > 3. As we 
move to higher 6-jet multiplicity, the signal distribution begins to stand out clearly from BG, 
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Figure 6: Distribution of number of jets in the FP5 case study with mo = 3050 GeV, rrii/2 = 400 
GeV, A = 0, tan/3 = 30, fi > and m t = 175 GeV (where nig = 1076 GeV), versus various 
SM backgrounds. 
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Figure 8: Distribution of number of isolated leptons for the FP5 case study with m = 3050 
GeV, mi/2 = 400 GeV, A = 0,tan/3 = 30, /i > and m t = 175 GeV (where m g = 1076 GeV), 
versus various SM backgrounds. 



which is dominated at this point by ti production. 

Alternatively, we can move to higher jet multiplicity. In Fig. [TOj we again plot At but 
this time for n(jets) > 7 and n(b — jets) > 2. The signal emerges from the BG clearly above 
At ~ 1300 — 1400 GeV, and has an advantage over Fig. [9]d) in that a somewhat larger signal 
rate remains after cuts. For the case shown, by imposing At > 1400 GeV, we are left with a 
signal cross section for case FP5 of 11.1 fb, while BG from ti production is at the 1.5 fb level 
with a tiny contribution from QCD multi-jet production. In addition, the remaining signal is 
98% from gluino pair production, so is almost entirely from the hard component of the signal. 

5 Signal, background and sparticle mass extraction 

We will adopt the cuts of Sec. H] as our cut set C2: 
C2 Cuts: 

apply cut set CI 

n(jets) > 7 
n(b - jets) > 2 
A T > 1400 GeV. 

These cuts have been optimized for m-g ~ 1 TeV. Next, we plot in Fig. [H]the event rate after 
C2 versus m g along a line of FP region with fl^ h 2 ~ 0.11, with tan/5 = 30, A = and 
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Figure 9: Distribution in At (defined in Eq. 8) in events with n(jets) > 6 with varying number 
of 6-tags, for the FP5 case study with m = 3050 GeV, mi/ 2 = 400 GeV, A = 0,tan/? = 30, 
/i > and m t = 175 GeV (where nig = 1076 GeV), versus various SM backgrounds. 
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Figure 10: Distribution of At in events with > 7 jets and > 2 6-tags, from the FP5 case study 
with m = 3050 GeV, m 1/2 = 400 GeV, A = 0,tan/3 = 30, // > and m t = 175 GeV (where 



m g = 1076 GeV), versus various SM backgrounds. 



/i > 0. For m g ~ 700 GeV, < 103.5 GeV, so the region is excluded by LEP2 chargino pair 
searches. The solid blue curve denotes the signal rate after cuts C2, while the brown dot-dashed 
curve denotes SM BG. Signal rates are typically in the multi-fb regime, and exceed BG out to 
m g ~ 1500 GeV. 

Since the signal in Fig. [11] comes from nearly pure gg production, the total rate can be used 
as an absolute measure of the gluino mass. There are of course a variety of theoretical uncer- 
tainties which arise. One comes from how well-known is the absolute gluino pair production 
cross section. The value of &{gg) has been computed to NLO in QCD in Ref. [36], where it is 
shown that a variation in renormalization/factorization scale leads to an uncertainty in u{gg) of 
±11%. A further uncertainty arises from variations in the squark mass. Here, we are assuming 
decoupled scalars, so variation due to changes in rriq are expected to be small. Nonethless, we 
find that by varying rriq : 2 — 5 TeV, the cross section still varies by ±10%. Folding the NLO 
uncertainty in quadrature with the rriq uncertainty, we estimate the cross section uncertainty 
at ±15%, and plot the expected theory cross section variation as the blue dashed lines. 

At this point, it can be asked how well will we know the gluino branching fractions, upon 
which the signal rate also depends. Here, we remark that in the region with decoupled scalars, 
we are relying on a value of \x that is just right so that the neutralino LSP saturates the CDM 
relic density measurement. Small variations in \i about this region are found to lead to only 
small changes in the gluino branching fractions. This is shown in Fig. [T2J where we plot in frame 
a) variations in Q,^h 2 versus /i, and in frame b) variations in the dominant gluino branching 
fractions. In the plot, we adopt as usual the case study FP5, and vary fi by adopting the 
non-universal Higgs soft mass model[10] in Isajet, which allows use of mSUGRA parameters, 
but also independent variation in the fi and rriA parameters (we keep fixed). 
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Figure 11: Cross section after CI plus > 7 jets, > 2 6-tags and At > 1400 GeV, for various 
points along the HB/FP region with ^l^Ji 2 ~ 0.11 with Aq = 0,tan/? = 30, \i > and m t = 175 
GeV, versus rrig. We also show a band of the theoretically expected uncertainty of our results 
due to variations in factorization/renormalization scale and variations in rriq ~ 2 — 5 TeV. We 
also show the level of expected SM background. 



It might also be argued that the event rate depends on the value of tan/3 that we have 
selected for our case study. In fact, since scalar masses have decoupled, b and r Yukawa 
coupling effects are tiny, and the variation of the signal after cuts C2 with tan (5 is comparatively 
negligible, as long as we require that the /i value be fixed so that one obtains the relic density 
^Zi^ 2 ~ 0.11. This is illustrated in Table HI where we plot signal rate after cuts C2 for 
tan/3 = 10, 20, 30, 40 and 50. In each case, the value of m 1( / 2 is fixed at 400 GeV, but mo is 
chosen so that the correct relic density is obtained. The resulting cross section after cuts C2 
shows only a ±6% variability. Meanwhile, variations in the Aq parameter again mainly affect 
the scalar sector, but since these decouple, the effects should again be small. 

A further consideration is to ask how well we really know our background estimates. At 
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Table 4: Cross section after cuts C2 for HB/FP cases with mi/ 2 = 400 GeV, A = 0, /i > and 
m t = 175 GeV. We list the m value required to give h 2 ~ 0.11 for different tan/3 values. 
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Figure 12: a) Variation in neutralino relic density with variation in \i for case study FP5. In 
b), we show variation in dominat gluino branching fractions versus ji. 



1(3 



this stage, the answer is difficult to know, and depends on several factors, including how well 
the selected event generator, Isajet, models SM backgrounds. If indeed ti production is the 
dominant BG, then the plethora of ti events produced at the LHC will allow detailed study of 
this reaction, so that the distributions will be well-known from data. Better theory modeling- 
such as inclusion of exact matrix elements for extra jet radiation[H]- will also help. Likewise, 
it can be expected that W + jets, Z + jets and QCD backgrounds will also be well-studied, and 
the high n(jet) and high At tails will be better known due to actual collider measurements. 
In any case, we try to make a rough estimate by simply assuming that our event generator 
background is known to ±100%. We add and subtract this BG uncertainty to our theory 
curves in Fig. [II], with the resultant band being denoted by orange dashed lines. 

At this point, we can try to estimate the precision with which the gluino mass can be 
extracted from total cross section measurements. We show in Fig. (TT] as data points the error 
bars expected from measuring the total cross section after cuts C2 with an assumed 100 fb _1 
of integrated luminosity (red data points). A simple estimate of the uncertainty can be gained 
from the intersection of the upper and lower limits on the statistical cross section measurement 
with the band of theory uncertainty. Using this method, we find that points 2-6 yield a gluino 
mass measured in the range of ±8%, as shown in Table [TJ The precision will increase or decrease 
depending on the ultimate uncertainty ascribed to the BG by the experimental groups. It would 
also decrease if an NNLO computation of gluino pair production is made. Note that even if 
the statistical error bars drop to zero (infinite integrated luminosity), the theory uncertainty 
still gives ~ 7% uncertainty. FP1- which is below the LEP2 excluded boundary- is difficult 
to measure because the projected theory curves level off for lower values of m g . This is just a 
result of the fact that we optimized cuts in the 1 TeV m g region. A better optimization with 
softer cuts would need to be performed to extract these lower gluino masses. For m g ~ 1300 
GeV, another optimization would be needed with harder cuts. Here, the absolute gluino pair 
event rate is dropping, so we expect a rate-based measurement of m g would be more challenging 
and perhaps not feasible in this higher mass region. 

6 Leptonic signatures 

While the analysis in Sec. \5\ focussed on lepton-inclusive signals, it is also useful to make use 
of the isolated lepton content of the signal. We expect events containing multiple isolated lep- 
tons to have somewhat reduced jet multiplicity compared to events with zero or one isolated 
lepton. To proceed with the multi-lepton channels, we retained cuts CI and examined the At 
distribution maintaining n(b — jets) > 2 but requiring n(jets) > 4 or 5. The distribution in 
At for n(jets) > 4 is shown in Fig. [T3J. Here we see signal emerging from BG for At > 1200 
GeV. (The plot using n(jets) > 5 is similar, but with lower signal and BG rates.) Hence, we 
adopt cut set C3 for events with 2 or more isolated leptons: 

C3 Cuts: 

cuts set CI 
n{isol. leptons) > 2 
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Figure 13: Distribution in At of events with cuts CI, n(leps) > 2, n(b — jets) > 2 and 
n(jets) > 4 for the FP5 case study with m g = 1076 GeV. 



In Fig. dH we show the signal rate of various multi-lepton topologies versus m g for FP 
cases with Aq = 0, tan/3 = 30 and /i > 0. The zero and one lepton topologies use cuts C2, 
while the same-sign (SS) dilepton, opposite sign dilepton (OS) and trilepton rates use cuts C3. 
We see that there should be consistent signals above SM backgrounds in all the various multi- 
lepton channels for much of the mass range of m g . Same sign lepton events would establish the 
Majorana nature of the gluino [42] . 

It is also well-known that kinematic information on neutralino mass differences can be 
gleaned by examining the invariant mass distribution of opposite-sign/same flavor dilepton 
pairs (OS/SF)[43]. We plot in Fig. [15] the invariant mass distribution for case study FP4. The 
HB/FP region is characterized by the fact that — < Mz and by m^r — < Mz, 

so that two-body spoiler decays of Z 2 and Z 3 are closed. We then expect two mass edges in 
the m(£ + £~) distribution: in the case of FP4, one is at — = 53.8 GeV and another at 
m^ 3 — = 75.1 GeV. Indeed, the double mass edge structure is becoming visible in the il% 
distribution with 100 fb" 1 of data as shown in Fig. [15] 



7 Conclusions 

In this paper, we have examined the sort of collider events to be expected at the CERN LHC 
for SUSY models in the HB/FP region of the mSUGRA model. We found that by requiring 
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Figure 14: Event rates for zero and one isolated lepton events after cuts C2, and OS and SS 
dileptons and trileptons after cuts C3, versus m- g . Zero and one lepton backgrounds are shown 
with the same color as the signal. We found no backgrounds events to OS, SS and 3/ to the 
cross section level shown. 




Figure 15: Distribution in OS/SF dilepton invariant mass for case FP4 using cuts C3. Two 
mass edges are becoming apparent at a luminosity of 100 fb _1 , in addition to the Z peak. 
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high jet and 6-jet multiplicity, and a high effective mass cut, a rather pure signal emerged 
from a dominantly tt SM background. Since the signal came almost entirely from gluino 
pair production, and the decay branching fractions were fixed by assuming the neutralino relic 
density saturated the WMAP Q^h 2 measurement, the total signal rate could be used to extract 
an estimate of the gluino mass. Factoring in theory uncertainty on the total cross section and 
a ±100% error estimate on remaining background, we found that rrtg could be measured to a 
precision of about 8% for 100 fb _1 of integrated luminosity. This was our central result. 

We note here that our conclusions apply more generally than to just the HB/FP region of 
the mSUGRA model. The key assumptions needed for our analysis are that: 

1. The flavor/CP conserving MSSM is the correct effective theory of nature at the weak 
scale, with the lightest neutralino as LSP. 

2. We assume gaugino mass unification, as occurs in many SUSY GUT and string models. 

3. We assume that scalars- the squarks and sleptons- decouple due to mass values beyond 
the few TeV level. This leaves just the various gluinos, charginos and neutralinos con- 
tributing to LHC collider events. 

4. The value of /i is fixed by the requirement that the relic abundance of Z\ saturates the 
WMAP measured value. This, along with gaugino mass unification, fixes the sparticle 
branching fractions to their assumed values. 

If these conditions are fulfilled, then the methods presented here should allow for a gluino mass 
extraction if m- g is in the mass range of ~ 700 — 1300 GeV. We note here that our result 
depends on the sparticle branching fractions being fixed to values near our calculated results. 
These in turn depend on the above assumptions being fulfilled. Thus, our study should be 
applicable to other heavy scalar situations, not only the HB/FP region of mSUGRA. Recently, 
such models have received renewed attention in light of FCNC constraints, and some string 
theory motivated models have produced heavier scalar spectra[H]. Our considerations also 
apply to the low scalar mass regime of split SUSY models {13], where the gluino decays promply 
inside collider detectors. 

In addition, we note that the signal can be separated as to its isolated lepton content. 
Typically, for each additional isolated lepton, there should be on average 1.5 less jets per event. 
The OS/SF dilepton mass distribution embedded in the hard signal component should exhibit 
mass edges at m% — and also at — m-x , which are distinctive of this scenario in which 
the LSP is a mixed bino-higgsino particle. The same mass edges should appear in the clean 
trilepton channel originating mainly from chargino- neutralino production (the soft component), 
as shown in Ref. [30]. The mass-difference edges, along with the absolute gluino mass, may 
provide enough information to constrain the absolute chargino and neutralino masses (including 
the LSP mass), under the assumptions listed above. 
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